Theoretical analysis of spin qubit decoherence near optimal working points 
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We investigate suppression of various decoherence mechanisms at critical points variously termed 
optimal working points (OWPs) or clock transitions (CTs) of electron paramagnetic resonance 
(EPR) transitions in Si-based spin qubits, where the frequency - field gradient is zero (df /dB = 0). 
New experiments by Wolfowicz et al. arXiv:1301.6567 (2013)] have demonstrated suppression of 
both instantaneous diffusion, as well as spin diffusion at these points, yielding orders of magnitude 
enhancement in coherence times: up to seconds, without extrapolation to zero-duration pulses as 
was until now necessary. We obtain simple analytical expressions and compare with full numerics, 
explaining the form of decoherence suppression for bismuth donors in Si (Si:Bi). It is demonstrated 
that behavior at the OWPs cannot be fully analyzed in terms of the behavior of df/dB, especially 
in the case of spin diffusion, indicating the need for an understanding at the level of the underlying 
microscopic processes entangling the central spin system and the spin bath. 

PACS numbers: 76.30.-v, 76.60.Lz, 03.65.Yz, 03.67.Lx 



i 

G 

O 

o 



> 
o 

CN 
O 
CO 



I. INTRODUCTION 

Coupled electronic-nuclear qubits have occupied a 
prominent role in the quest for robust and scalable quan- 
tum computing ever since the seminal proposal of Kane-i 
for exploiting the long-lived coherence of nuclear qubits, 
in conjunction with the possibility for fast manipulation 
of the electron spin. Most work has considered phospho- 
rus ( 31 P) donors in Si (Si:P)ii~— In recent years, there 
has been growing interest in 209 Bi, the deepest Group V 
donor in Si (Si:Bi)ji2r— Various approaches for controlling 
decoherence have attracted considerable effort: dynami- 
cal decoupling techniques j2Sr^£ the use of so-called "clock 
transitions" (points at which the frequency-magnetic field 
gradient df/dB vanishes) in systems as diverse as flux 
qubits, 30 nuclear spin qubits^Ir— and electronic spin qubits 
of Si:Bi; 17,20 finally, studies of isotopically pure samples 
have indicated that electronic Ti coherence times can be 
of order of seconds 34 while for nuclear spin qubits T2 can 
reach minutest 

Isotopic purification is effective in reducing decoherence 
because for donors in Si and other related systems, deco- 
herence is primarily due to host impurities - 29 Si in the 
case of Si, since the 29 Si isotope (with natural abundance 
p = 4.67 %) has nuclear spin 1 = 1/2 and introduces spin 
diffusion, mainly by spin flip-flops. The availability of sam- 
ples comprising almost pure 28 Si (which has no spin) might 
appear to offer a definitive solution by eliminating the de- 
cohering 29 Si spin-bath entirely. However, the previous 
measurements of intrinsic electronic T2 times of order sec- 
onds in purified samples were extrapolated from the use 
of electron paramagnetic resonance (EPR) pulses of zero 
duration, in order to suppress instantaneous diffusion, a 
result of donor-donor interactions which otherwise limits 
T2 to tens of milliseconds. Unfortunately, pulses of zero 
duration preclude a reasonable implementation of qubit ro- 
tations (e.g. a Hahn echo sequence) and restrict the possi- 
bilities for pulse sequences relevant to quantum information 
processing. 

Now, a new experimental study of Si:Bi at CTs 3 ^ demon- 



strated coherence times on macroscopic timescales with- 
out having to perform partial rotations and even for full- 
duration 7r-pulses. In addition, in natural Si, a factor of 
100 enhancement in T2 was observed at an OWP, to about 
93 milliseconds. It has been previously demonstrated for 
Si:P that 29 Si nuclei suppress donor-donor flip-flops^ As 
pointed out in 3 ^, the optimal situation might arise from a 
combination of partial enrichment and CTs. 

Although CTs or optimal working points (OWPs) are 
of increasing interest, understanding of the underlying 
physics for spin decoherence is still lacking. Theoretical 
analysis is usually framed in terms of the qubit's reduced 
sensitivity to externally induced and uncorrelated magnetic 
field noise at the df/dB = points; the T2 behavior is then 
related to the curvature of the f(B) function 31-33 including 
higher order behavior such as d 2 f/dB 2 . In previous studies 
using nuclear spin qubitsj 3 ^— these points are termed Ze- 
FOZ (Zero First Order Zeeman Shift) points, rather than 
CTs, while in the earliest studies on flux qubits^ - the term 
OWPs is used. For convenience, we refer to them as OWPs. 

Parallel to these studies, the cluster correlation expan- 
sion (CCE) has provided a versatile and powerful frame- 
work for realistic simulation of the decoherence of a cen- 
tral spin in a spin bath42r— . For the important problem 
of spin diffusion, the key physics leading to decoherence 
is entanglement between the bath and central spins i 43 i 44 
Thus, analysis of decoherence in terms of sensitivity to a 
randomly fluctuating external magnetic field is at best an 
emergent behavior. We show in Fig. [T] by comparison with 
full CCE numerics, that the suppression of spin diffusion 
due to indirect flip-flops by bath spins cannot be described 
by df/dB, except locally. 

In the present work, we obtain simple analytical expres- 
sions which clarify the physics of OWPs and explain the 
main experimental features. They are derived from a sim- 
plification of the CCE formalism (which we term "uCCE" 
for unperturbed CCE) using small clusters. 

The analysis suggests that spin diffusion at OWPs is 
dominated by a relatively small proportion of strongly cou- 
pled spin clusters (i.e. for which the coupling to the central 
spin system is much stronger than the intrabath dipolar 
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FIG. 1: (color online) Comparison between numerical T2 times 
due to spin diffusion, calculated using the cluster correlation 
expansion and from df/dB. It can be seen that the behavior 
around optimal working points (OWPs) cannot be described by 
df/dB, except locally: the df/dB curve can either be scaled to 
agree with the asymptotic behavior (blue) or the inner region of 
the OWP (red) but not both. Thus, analysis of the underlying 
microscopic processes is essential to yield the full magnetic field 
dependence of 22(B). 
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FIG. 2: (color online) Electron spin coherence times T2, due to 
spin diffusion, at five Si:Bi EPR optimal working points (OWPs) 
evaluated for a Hahn echo sequence, FID and with our analyti- 
cal formulae. Note that the analytical formula reproduces accu- 
rately the different shape of OWPs 1-4 relative to OWP 5. The 
latter is an NMR transition at high magnetic fields, but is an 
EPR transition at the OWP 5 point of B = 0.37 T. The green 
line indicates the result of convolving the behavior with a Gaus- 
sian of FWHM 2.1 G to allow for inhomogeneous broadening 
in natural Si. This reduces the maximum T2 to about 100 ms. 
The intrinsic T2 maximum for a single bismuth in natural Si 
at OWPs 1-4 is of order seconds, but experiments sample an 
ensemble of bismuth donors, with a 0.42 mT spread in field. 



coupling between the flip-flopping bath nuclei). We obtain 
the shape of the OWP "spikes" {T 2 (B) plots) which for a 
single central spin system in a bath are limited to a maxi- 
mum T 2 ~ 10 s. However, the upper limit to T 2 at OWPs 
for an ensemble of donors is of order 100 msj^ which here 
we attribute to inhomogeneous broadening. The uCCE ex- 
pressions make clear the nature of the observed strong sup- 
pression of instantaneous diffusion and other donor-donor 
processes at OWPs in Ref. [H. We conclude that while the 
T 2 (B) due to indirect flip-flopping cannot be well described 
over the full range by a simple df/dB curve, in contrast, in- 
stantaneous diffusion, which is related to a simple energy 
shift, does follow an approximate (df/dB) 2 behavior, as 
seen exp er iment ally 2 1 1 36 This study does not eliminate the 
need for careful, high-order, large-cluster CCE numerics, 
which yield highly accurate quantitative results. 48 How- 
ever, a key feature of the CTs is that the T 2 (B) curves 
can vary by over three orders of magnitude over an ex- 
tremely narrow magnetic field range (tens of G) as seen in 
Fig. [2] Thus, we aim primarily to understand the underly- 
ing microscopic processes and to estimate T 2 (B) to within 
a factor of 2 or so. 



This paper is organized as follows. In Section II, we 
develop a model of spin decoherence derived from the CCE 
and use it to obtain analytical expressions describing the 
behavior of OWPs. We apply the theory to nuclear spin 
diffusion in a bath of spin- 1/2 nuclei, and further consider 
Si:Bi as the central spin system. In Section III, we present 
an analysis of donor-donor decoherence processes. 



II. MICROSCOPIC THEORY OF NUCLEAR SPIN 
DIFFUSION AT OPTIMAL WORKING POINTS 

A. Cluster correlation expansion 

The methodology of the CCE is extensively detailed else- 
where both for a single central spin ; 43 ' 44 as well as an 
ensemble4^ In calculating successive terms in the expan- 
sion, the Hamiltonian including a central spin system is 
diagonalized for sets or "clusters" of bath spins of varying 
sizes until convergence is achieved. The spin echo or free 
induction decay (FID) intensity computed for the largest 
cluster (containing the entire bath) as a function of time, 
t is denoted by C(t). It corresponds to the off-diagonal el- 
ement of the central system's reduced density matrix, nor- 
malized to unity at t = 0. The expansion takes the form 
of 

C{t)=Y[C K (t). (1) 

K 

where JC denotes a cluster, and the product is over all clus- 
ters in the bath. Using Eq. ([T]), and defining £jc(t) as the 
echo calculated for spins contained only in K,, a recursive 
expression for each cluster term (Cic{t)) is obtained: 

C K {t) = C K (t)/ H £ K ,{t). (2) 

The exact echo C(i) is recovered in the limit of includ- 
ing all clusters in the expansion. However, for spin diffu- 
sion, convergence is obtained by including only those clus- 
ters with at most a pair of bath spins. This is due to 
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the large energy mismatch between the interaction among 
bath spins and the coupling of the central spin system to 
the bath4i A modified CCE was recently developed for the 
case when this energy mismatch is small.— More generally, 
the ^-cluster CCE is defined as 

£ (i) (t)= n c K (t), (3) 

\K\<1 

which includes clusters with at most I spins. We confirmed 
the expected convergence for spin diffusion as I = 2 — > 3. 

B. Nuclear spin diffusion 

Within the CCE framework, we consider a central spin 
system, which consists of an electron spin coupled to a 
nuclear spin, and interacting with a bath of other spins. 
The total Hamiltonian is written as 

-Htot = Ho + -ffint + -ffbath- (4) 

The first term denotes the Hamiltonian for the central spin 
system. In nuclear spin diffusion, the central spin system 
interacts with clusters of nuclear spins: 

^int = ^^SJ„i n , (5) 

n 

where J„ for the n-th spin is in general of tensor form 
(for anisotropic couplings). The final term represents the 
bath Hamiltonian including the Zeeman interaction with 
the external magnetic field B: 

(6) 

n n<rn 

where ujq — fj,B/h is the electronic Larmor frequency 
(p = 1.857 x 10 -23 JT _1 ), Sn is the ratio of the nuclear to 
electronic Larmor frequency and r nm denotes the relative 
position vector of bath spins at lattice sites n and m. The 
components of the dipolar tensor D (r) are given by 

Mo ^/i 2 / 3r«rA 

where <5y denotes the Kronecker delta, /io = 4ir x 
10~ 7 NA~ 2 and i, j — x, y, z. 

C. uCCE approximation 

While for the full CCE we would diagonalize H to t, here 
we introduce a simplified model. We first diagonalize Hq: 

H \i)=E^\i). (8) 

Typical values of Eq are in the GHz range for electronic 

spin states. We then diagonalize h^' — (i\ (-ffint+ffbath) 
in other words the bath/interaction Hamiltonian averaged 
over the state of the central spin system. For nuclear 
spin diffusion, this is a reasonable approximation given 



Hq ^> -ffint,-ffbath as the bath and interaction terms are 
on the kHz-MHz scale. Previous studies have employed 
the simplification of a bath evolution conditional on the 
state of the spin; 46 i 47 for systems where the central spin 
state is a simple electronic Zeeman state, whereas here we 
extend this to a more complex system comprising a highly 
mixed nuclear-electronic state, in the OWP regime. The 
key idea is still that the central spin regulates the bath 
evolution differently, depending on whether it is in the up- 
per or lower states. However, here it is a strongly mixed 
state of 'up' and 'down' electronic/nuclear states, but the 
GHz scale of the central spin eigenfrequencies means that 
the bath 'sees' only the fast-averaged effective magnetic 
field. We proceed with diagonalizing the Hamiltonian ftW 
(and it is convenient here to use as a basis the states \k B ), 
which are eigenstates of the non-interacting bath). The 
eigenstates of hS % > are: 

h®\k i )=E®\k i ) (9) 

One immediate effect of the approximation is to in effect 
reduce the interaction with the bath to Ising form H mt = 
J n S z I* in Eq. ©. Since here we neglect coupling between 
the nuclear spin of the central spin system and the bath, 
in , we simply replace the operator S z by the electronic 
polarization for each eigenstate (the difference between the 
probability for an up and a down electronic spin for the 
central spin system): 

Pi (w ) = 2 (»| S z \i) . (10) 

Near an OWP, Pi (luq) varies rapidly and by orders of mag- 
nitude with the external magnetic field B. For an unmixed 
state, Pi (wo) = ±1. For an arbitrary \u) — > |Z) EPR transi- 
tion, \P U — Pi \ ~ \df /dB\. Nevertheless, it is worth noting 
that df/dB is a property of Hq (the central spin system) 
while P u and Pi refer to the interaction term few, As the 
former is of predominantly Ising form, dependence on P u 
and Pi can locally be related to df/dB. 

We now consider a Hahn echo sequence of an arbitrary 
\u) — > \l) EPR transition. At time t — 0, the system is 
prepared by means of a ^-pulse in the state \t = 0): 

|t = 0) = ^(|u) + |i))|fc B ), (11) 

the bath in initial state |fcs). A 7r-pulse is applied at t = 
r, then after further evolution for a second period r, the 
corresponding Hahn echo for the /C-th spin cluster at t — 2r 
is given by 

mzL (£> E c kfj,r e~ u£ ^' , 

k,k',j,j'=l,\K\ 

(12) 

where 

C u k ^ hf {k B ) = (j'Jji){k u \ji){fM{k u \k B ){k B \k[)^ 
£J$ l M> = E)-E}+EZ-E l k ,. (13) 

Note that if we calculate J2k B ^Hahn ^B,t), we can use 
the completeness of the basis of non-interacting bath states 
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to simplify Eq. ([12]) . The intrinsic FID is simpler: 



2x2 Hamiltonian analytically, obtained from 



r u 

*-F 



FID (^J k B ,t) 



2fc,fe'=l,|/C| 

<WC/;IWIIWbI*») 



(14) 



We note that the above expressions are valid for any cluster 
size, |/C|, not just two-clusters. The key assumption above 
(in addition to neglecting the non-Ising part of the interac- 
tion between the central spin system and the bath) is that 
since Po 3> Pbath we also neglect the "back-action" of ftW 
on the central spin states: P, (u>o) is evaluated without con- 
sidering the interaction term. We refer to Equations (Tl"2")) 
and (0 as uCCE ("unperturbed-state CCE") simply to 
distinguish them from the numerical results of the CCE 
arising from full diagonalization of Ptot ■ Below we investi- 
gate primarily two-cluster uCCE (i.e. \JC\ =2) but test full 
numerics with \JC\ — 3. One advantage of the uCCE is the 
significant computational saving. For bismuth donors and 
similar systems in 29 Si baths, our comparisons suggest they 
are extraordinarily accurate to within 0.1 G of the OWP 
point. In general, they follow in detail the equivalent be- 
havior obtained from full diagonalization which includes 
the non-Ising part of the system-bath interaction. 

A potential advantage for an OWP investigation, is that 
the the C^ n {JC,k B ,t) in Equations and (fTU) typi- 
cally oscillate on the kHz scale, while the comparable CCE 
quantity oscillates on the GHz scale. Near OWPs, we may 
wish to propagate C(t) on a millisecond timescale at one 
magnetic field, while at a nearby field we may need to prop- 
agate on the timescale of nearly a minute; i.e. while E% 
and E l k , are individually on GHz timescales, the E k — E l k , 
differences in Equations (|12j) and (|14j) are comparatively 
small. The elimination of this potential numerical prob- 
lem may make propagation of the uCCE on macroscopic 
timescales more stable. 



D. Spin diffusion from spin-1/2 flip-flops 

In natural Si, the limiting decoherence mechanism at 
low temperatures is spin diffusion in a bath of spin-1/2 
29 Si impurities. Experimentally, decay of the Hahn echo 
for electronic spin qubit systems is typically fitted to 
exp[— t/T 2 — (t/TsT>) n ], where Tsd < 1 ms characterizes 
the spin diffusion with n ~ 2 ; — and other processes rep- 
resented by T 2 . Although T 2 3> Tsd, for convenience, we 
refer to our calculated decay time-scale as T 2 below, noting 
that we are in fact investigating spin diffusion. 

Conveniently, Eq. (fl"4|) reduces to simple algebraic ex- 
pressions: there is no need at all to undertake any numer- 
ical diagonalization in order to compute T 2 . This suffices 
to reproduce most of the complexities of the Si:Bi OWP 
behavior described below. 

For spin 1/2 systems like 29 Si, bath states \ks) corre- 
spond to |t|), ||t), |ft) and We note that |ft) and 
|H) cannot flip-flop, thus in our simple (Ising-only system- 
bath interaction) model, provide no decoherence. For the 
flip-flop states |f|) and ||t), we just diagonalize a simple 



J=l,2 

. (15) 

corresponding to a single cluster of two 29 Si spins inter- 
acting with dipolar coupling strength Ci 2 , and each cou- 
pled to the central spin system with strength J\ and J 2 
respectively. Let i denote either u or I (ie the upper/lower 
states of the EPR transition). We introduce the parame- 
ters = (Jf - JK)P and e t = ^{{Cf 2 f + (Cf ) 2 ). For 
simplicity, below we drop the cluster label /C, but it is un- 
derstood that all quantities are calculated for a single 29 Si 
spin-pair, of which there are N ~ 10 4 in the bath. The 
eigenvalues are AE± = — C\ 2 ± a. Defining Q. — Ci ± 
the eigenstates are 



4 lit) |i) + bf |U) \i) 



(16) 



with of = C 12 [(C ± ) 2 + C\ 2 Y^ and bf = C±[(C±) 2 + 
C 2 2 ] _1 ^ 2 . The eigenstate overlaps (e.g. (k u \k[)) in Equa- 
tions (fl"4]) and (fT2j) are trivially evaluated from the af and 
bf. We thus rewrite: 

2£« ID i (/C,tl,t) - D+e-«<~- c * + D- t e + *«'- c * 

+R+e-^+^ t + R- l e +l ^+^ t (17) 



where 



P 



[(C? 2 + (g)*)(C& + (&)*)] 



[{a 



(Q) 2 )(c 2 12 



m 2 ) 



(18) 



The 



for 



expression ^ ^ Hahn 
gous, but much more complex 



is analo- 
In addition, 



£Fro(£>tt)*) = £fid(£,4-M) = !■ Thus we can replace 
the average over bath states by the simple approximation 

The suppression of decoherence at OWPs is considerably 
clarified by Eq. (|T4"|) . It gives the decay for a single cluster 
(the decay for the whole bath is CpJ^ (t) = Yl K £piD ) 
and comprises four oscillating terms. At the OWP, e u = q 
so the terms dependent on D ± cease to oscillate and stop 
contributing to decoherence. The P terms continue to os- 
cillate, rapidly. However, with some manipulation it is easy 
to show that at the OWPs, while D ± -t 1, the P ± 0, in 
other words, the amplitude of the rapidly oscillating terms 
vanishes. The decay behavior close to OWPs is entirely 
dominated by the P terms. 



E. Application to bismuth donors in Si 

The Si:Bi system has an electron spin S =1/2 coupled to 
a large nuclear spin / = 9/2 by a strong hyperfine coupling 
constant, 4~ = 1475.4 MHz, 



Po — ^0 



(s z + s Bi i z ) 



AI.S, 



(19) 
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where 5bi = 2.486 x 1CP 4 is the ratio of the bis- 
muth nuclear to electronic Zeeman frequencies. Measure- 
ments suggest isotropic couplings dominate the system- 
bath interaction. 2 — The eigenstates \i) correspond to i = 
1,2... 20. The quantum numbers i are used to identify 
EPR transitions as shown in Fig. [5J thus OWP 2 corre- 
sponds to the experimentally investigated CT in Ref. 
(| 14) — > 1 7) in the above notation). 

The only parameter required to calculate a T 2 time, in 
addition to a set of coupling strengths G\ 2 , J\ and J 2 , is 
the polarization Pi of the central spin eigenstates. The po- 
larization is strongly field dependent near OWPs. The en- 
ergy eigenstates \i) of Si:Bi near OWPs are strongly mixed 
and are not well approximated by the usual Zeeman ba- 
sis \ms,mi). The mixed eigenstates may be labeled in 
an adiabatic basis and correspond to doublets of constant 

m = 777 5 + 777/: 

|t) = \±,m) =a± \±\,mT\) + Pt\T\,m±\) ■ (20) 

All quantities are easily obtained analytically^ The elec- 
tron polarization as a function of magnetic field is 

Pi m = n m {n 2 m + 25 - m 2 y 1/2 , (21) 

and Cl m = 777 + ^ (1 + (5bi) ; m is an integer —5 < m < 5. 
For the experimentally observed transition in Ref. l36l . m — 
— 1 for the upper state, while m = —2 for the lower state. 

The simple expression Eq. (flTT) matches accurately to 
(within ~ 1% or so) the decay behavior of a correspond- 
ing numerical CCE calculation, whether far from OWPs 
or within the ~ 100 G width of OWPs 1-4. The only dis- 
agreement is well within 0.1 G of the OWP, where the full 
numerics suppress the divergence in T 2 . 

Figure |3] illustrate our results. We plot only uCCE FID 
results from Eq. (jTTf since the results from full CCE FID 
numerics were completely indistinguishable in the plots dis- 
played. The only exception is the regime within 0.1 G of 
the OWP: the uCCE results yield T 2 — > 00 at the actual 
OWPs, whereas the full numerics saturate at T 2 ~ 10 s. 
This is beyond the range shown in our main result in Fig.[5J 
For natural Si, this discrepancy does not affect the mea- 
sured results, since inhomogeneous broadening of 4.2 G 
means that the weight of the region within 0.1 G of the 
center is rather small. 

Failure of convergence indicating the need for inclusion 
of higher-order corrections (larger clusters) manifests as 
oscillatory (non-decaying) behavior in our calculated Halm 
echos as a function of time, at longer time scales. This does 
occur as we approach an OWP. A simple expedient (for a 
non-converged calculation) is to consider the decay rate at 
short times as it can provide a reasonable indication of the 
magnitude of T 2 : inclusion of higher-order effects prolongs 
the calculation rather than radically changing T%. How- 
ever, fortunately, for reasonable parameters, the calculated 
FID does not fail and always displays non-oscillatory decay 
at all timescales. Thus, within about 50 G of the OWP, 
72-Hahn begins to fail; it ceases to decay at long times, 
and exhibits oscillations heralding loss of convergence. We 
can extend convergence times by using larger clusters, but 
cannot prolong this to the seconds timescales relevant to 
OWPs due to computational costs. In the experiment, the 
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FIG. 3: (color online) Comparison between uCCE expressions 
(red) and CCE (black), showing excellent agreement up to 
within about 0.1 G of the OWP. We note that the width of 
the OWP 1-4 peaks is of order 100 G, as shown in Fig. [5J so 
uCCE fails over a very narrow region. Within about 0.01 G of 
the OWP, the uCCE T 2 time diverges while the CCE is limited 
to a maximum T2 of order 10 s. 



FID is affected by inhomogeneous broadening as it arises 
from an ensemble of bismuth donors, so the Hahn echo 
is essential. Our calculated FID in contrast is the intrin- 
sic decay of a single bismuth donor. This intrinsic FID 
(termed iFID), away from OWPs, differs negligibly from 
the Hahn echo and provides a reasonable estimate of the 
decay timescale throughout. Away from OWPs, T 2 _Hahn 
differs from T2_iFio by less than 10% at all timescales 
and there are no convergence problems for either. In 
Fig. [2] we show T2_Hahn values near OWPs obtained from 
the short time behavior of the Hahn echo, showing that 

— Hahn 

~ 2T2_iFiD, to an excellent approximation at all 
the OWPs. 



F. Approximate T2 expressions near OWPs 

We now derive analytical T 2 formulae, assuming 
T2_Hahn — 2T2_iFiD represents a reasonable approximation 
to the decay timescale. From Eq. p7|) . from the short time 
behavior we can infer that close to an OWP each cluster 
K contributes \C u ^ l (K.,t)\ w exp [-(t/T^) 2 ] where 




The above result is obtained by noting that near an OWP, 
&u —> ej (easily shown) and (D^+D^) —> 1, while R^ —> 0. 
The total T 2 is obtained by summing over all clusters: 



Equation (j2"3"]l predicts the intrinsic FID behavior. Dou- 
bling it represent an estimate of the Hahn echo decay. The 
above expression is displayed for OWPs 2, 3 and 5 in Fig. [2] 
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showing excellent agreement (violet dot-dashed lines) very 
near the OWPs. 

In Fig. 0] we use Eq. (l2"3")l to estimate the strength of 
the ZC-th 29 Si spin pair's individual contribution to deco- 
herence (its 1/(T^ C ) 2 ) as function of the hyperfine coupling 
to the central spin. Note that J\ and J 2 always appear 
as a difference in the equations, thus \J\ — J 2 \ is the rel- 
evant parameter. A striking pattern appears: the spins 
are grouped into lines of constant C12, corresponding to 
adjacent spins, second nearest neighbor and so forth. Sur- 
prisingly, for the stronger spins, the 1/(?2 C ) 2 is completely 
independent of the hyperfine coupling to the central spin. 
The origin of this behavior is clear from Eq. (|22p . For large 



I Ji - J2I > C12, the amplitudes R^i+R u i °c (j^Tj^ while 

the frequency of the oscillation (e u +ei) 2 cx (J\ — J2) 2 , elim- 
inating the dependence on the coupling. This conclusion 
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FIG. 4: Estimated weight of each spin in the bath (a total of 
10 4 flip-flopping pairs, as indicated by their individual contri- 
bution the total (I/T2) 2 , plotted as a function of the strength of 
the hyperfine coupling to the central spin system. We see that 
the 29 Si spins in the bath fall neatly into groups corresponding 
to different C12: here n = 1, 2, 3.. indicate nearest neighbors on 
the lattice, second-nearest and so forth. This behavior is valid 
only close (within ~ 100 G) of the OWPs. The scale of T 2 is 
set by the comparatively small N ~ 10 2 set of strongly-coupled 
spins in the red box. Data corresponds to B — 0.07978 T, i.e. 
about 0.5 G from the maximum of OWP 2. 



must be approached with some caution. Eq. ([22]) is valid 
only near the OWP and refers to short time decay behavior. 
Nevertheless, numerical test s do indicate insensitivity to 
the hyperfine coupling in the strongly coupled spins. The 
pattern persists unchanged for nearby B values. 

What is striking is that in Fig. [3J the N ~ 10 2 strongest 
spin pairs suffice to set the scale of T 2 . We have tested 
this by running a two-cluster CCE calculation with just 
120 spins. The FID yield T 2 ~ 140 ms with all - 10 4 spins 
and T2 ~ 160 ms with just these spins. All these spins 
correspond to | J\ — J 2 \ >• C12 and also to the strongest 
values of nearest-neighbor dipolar coupling C12 = 0.23 x 
10~ 2 MHz. 



We can use this to derive an approximate formula for 
decoherence near OWPs: 



T 2 x 10 d « = 



16 V|P u (w )fl(w )| 



C 12 Vn\Pu(uo)- Pi(uo)\' 



(24) 



The 10 3 scaling simply ensures an answer in milliseconds 
as all our energy units are in MHz. From the above, we es- 
ti mate C 12 y/N ~ 0.024 MHz. Over the narrow OWP 1-4, 
y/\P u (wo) Pi (wo)| is nearly constant, thus the field depen- 
dence reduces to 1/ (P u — Pi). Since (P u — Pi) cx df/dB 
there is an approximate (df/dB) 1 dependence close to 

the OWP. However, the additional y/\P u (wo) Pi (wo)| term 
in the numerator accounts for the much larger T2 in OWP 
5, without the need for scaling or fitting parameters, indi- 
cating that the (df/dB) -1 does not represent the full field 
dependence. 

From Eq. ([2~4"|) , we infer that the intrinsic T2 depends on 
1/y/p where p = 4.67 % the percentage of 29 Si in natural 
Si. Given that the inhomogeneous line width also varies as 
1/y/p, the combined dependence should be linear ~ I /p. 
It is reasonable to assume that the T 2 limit of 100 ms in 
natural Si will reach the seconds timescale for very modest 
enrichment factors of 100. Above this level, donor-donor 
effects are the limiting factors. 

All the uCCE expressions indicate a divergence at 
the OWP. However, this presupposes that all spin 
clusters can achieve the condition P u (ojq) = Pi (u>o) 
at the same external magnetic field ujq. In reality, 
Hbath does very slightly perturb the H eigenstates in 
Eq. {5J. The hyperfine coupling to the bath pro- 
vides a small correction to the former and by implica- 
tion to the polarization, P u (ujq) — > P u (cjo + SB) where 

SB ~ ^|af| 2 - |^'| 2 ) (Ji - J2) /4 V In addition, there is 
a correction from the non-Ising part of the interaction, 
which mixes the bismuth eigenstates. Full CCE numer- 
ics diverges from uCCE within 0.1 G of the OWP. Within 
~ 0.01 G the two-cluster CCE results indicate a maximum 
T2 of 10 s, suppressing the divergence. This is the true limit 
for a singl central spin system in a bath. However there 
are other sources of spread in the effective external field 
due to the fact that the experiment samples an ensemble 
of bismuth donors, subject to variations in the local fields. 
One is the local variation in A, the bismuth hyperfine con- 
stant, due to lattice strain^ With a 200 kHz line width 21 
this represents also a < 0.1 G spread. More drastic is the 
line width of natural Si, arising from the sea of single 29 Si 
spins (more numerous than the spin pairs responsible for 
spin diffusion) which cause a variation in B among spin 
qubits. This is of order 0.42 mT. Hence we convolve the 
OWP 2 behavior with a Gaussian to simulate this effect. 
We calculate the time decay 



D(t) 



1 



Of 



e (t/T2(B)f dB 



(25) 



where w = 0.21 mT and T2 (B) is our calculated X2 values. 
D(t) is found to give a non-Gaussian decay (stretched ex- 
ponential) as observed and reaches its e _1 value at 100 ms. 
Figure [2] also shows the result (green curve): the maximum 
T2 attainable OWP 1-4 in natural Si is of order 100 ms. 
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The close agreement with the experimental value of 93 ms 
is fortuitous: our analysis aims only for accuracy of a factor 
of order 2 or so, but this gives an indication of the expected 
experimental behavior. 

When using uCCE expressions rather than full numer- 
ics, we imposed a cut-off in the maximum Ti = 10 s very 
close to the OWP and do not permit higher values. How- 
ever, the convolution is not at all sensitive to the manner 
this intrinsic limit of 10 s from CCE numerics is imposed: 
the uCCE divergence is extremely narrow, thus the uCCE 
remains a reliable approach. 



III. DONOR-DONOR DECOHERENCE 

For the purest isotopically enriched samples, decoher- 
ence by 29 Si bath impurities is eliminated and the coher- 
ence times are instead limited by donor-donor spin inter- 
actions. The most destructive to previous experiments has 
been instantaneous diffusion. In a Hahn-like echo sequence, 
neighboring donor spins are also rotated by applied pulses: 
this is especially problematic for the second, 7r-pulse, as 
it leads to a temporally changing field whose effects are 
not removed by the echo sequencei 49 i 50 The other relevant 
processes are indirect flip-flops, whereby a central donor 
interacts with a bath of dipole-coupled donors, and direct 
flip-flops which involves the central donor flip-flopping with 
a bath donor. The latter is the limiting decoherence pro- 
cess. 



A. Instantaneous diffusion 

In this case, for uCCE, we revisit the key Hamiltonian 
in Eq. (|26[) which now has simply two donor spin systems, 
spins A and B (which can be hybrid electronic-nuclear 
states or more complex systems). 



fjA 



H, 



H, 



(26) 



where H m t — Dx.S A S A for coupling between donor-donor 
pair JC. We diagonalize H$ and Hq then consider the 
averaged interaction field (ia^bI H- lnt |za«b)- We consider 
the pulse sequence 



Ri {26 A )R*{2e B )-e 



iHi 



(27) 



meaning that the central spin A experiences a rotation of 
tip-angle B a (which in a Hahn sequence is 29 a = t) while 
the 'bath' spin B, unintentionally flipped by the pulse, un- 
dergoes a rotation by tip-angle 9 b (which in general need 
not be identical to 9a)- Here, R x represents a rotation 
about the a;-axis of the Bloch sphere. 

We assume the initial state to be -j= (|m) + \n)) A \m) B . 
Since both A and B have to be resonant with the pulse, 
their spin states must correspond to either of the states of 
the EPR transition i.e. m,n = u,l for \u) — > \l). The echo 



is given by: 



|£(2 7 



iT E m -E„- 



DP m (P n -P m ) 



sin 9 a + cos ( 



iTD(P m -P n ) 2 



sin 2 9b + cos 2 9b 



(28) 



where Dk = D for convenience. This equation is equiv- 
alent to Eq. (fH)) , but requires little further analysis. By 
inspection, we see that at the OWP, where P n = P m , the 
echo is suppressed quadratically in (P m — P n ). The spin- 
bath coupling strength D vanishes from the expression and 
in the second term, the tip-angle of the bath spin 9b also 
vanishes entirely, making the decay rate completely inde- 
pendent of the angle that the bath spin is (unintentionally) 
rotated through: 



\C(2r)\ 



1 



OWP 



COS( 



suit 



(29) 



For a 7r pulse, 29\ — tt and \£ (2r)| OWP = 1 for all B spins, 
regardless of coupling strength or 9b- 



B. Indirect flip-flops 

The case of indirect flip-flops is analogous to a pair of 
flip-flopping 29 Si's with the hyperfine couplings to nuclei 
replaced by dipolar couplings to a pair of donor electrons. 
The couplings to the central donor are scaled by the polar- 
izations P m and P n and the coupling between bath donors 
by (a m a n ) 2 . In addition, the eigenvalues of the interaction 
Hamiltonian arc shifted by +\c mn (1 — P m P n ). Like 29 Si 
spin diffusion, indirect flip-flops are suppressed linearly in 
(P m — P n ) at OWPs as discussed above. 



C. Direct flip-flops 

Direct flip-flops are not suppressed at OWPs and limit 
the electron Ti at low temperatures. The effective Hamil- 
tonian is H — Hq + Hq + Hint with 



Hu 



DS* A S% 



(S A S 1 



S A S BJ 



(30) 



For the initial state A= (|m) + \n)) A \m) B , the FID is given 
by 



\C(r)\ 



cos 



(31) 



and for a coherent superposition of \m) B and \n) B in the 
initial state, 



KM 



\ 



l + C os{i^ [(P m -P„) 



2a 2 n a 2 



(32) 



At an OWP, the (P m — P n ) 2 term vanishes, but this is in- 
sufficient to eliminate the oscillations in C (r) which cause 
decoherence. Thus, direct flip-flops are not suppressed at 
OWPs. 
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IV. CONCLUSIONS 

We present a theoretical analysis of spin decoherence 
near OWPs, which clarifies which decohering processes are 
suppressed and which are not. The method is a simpli- 
fication of the CCE approach, but applicable to general 
central spin systems which are not single spin- 1/2 systems 
(which have no OWPs). An advantage of the approach 
is that the central spin system need only be diagonalized 
once, thus saving significant computational time. In addi- 
tion, evolution for longer time periods is more numerically 
stable. Thus this simplification may find applications even 
in regimes without OWPs. Our expressions for T2 times 
are in agreement with full numerical diagonalization (ie ap- 
plication of the usual CCE method, including non-secular 
terms). 



In applying the theory to Si:Bi nuclear spin diffusion, 
we find that the behavior of T2 does not follow the form 
of df/dB except locally. Finally, we present an analysis of 
donor-donor decoherence effects which are known to limit 
T2 in isotopically pure samples. As demonstrated in recent 
experiments, we explain the suppression of instantaneous 
diffusion at OWPs and show that direct donor-donor 
flip-flops provide the limiting decoherence mechanism at 
low temperatures. 
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